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Observations of star motion, emissions from hot ionized gas, gravitational lensing and other 
tracers demonstrate that the dynamics of galaxies and galaxy clusters cannot be explained by the 
Newtonian potential produced by visible matter only [1-4]. The simplest resolution assumes that 
a significant fraction of matter in the Universe, dominating the dynamics of objects from dwarf 
galaxies to galaxy clusters, does not interact with electromagnetic radiation (hence the name dark 
matter). This elegant hypothesis poses, however, a major challenge to the highly successful Standard 
Model of particle physics, as it was realized that dark matter cannot be made of known elementary 
particles [4]. The quest for direct evidence of the presence of dark matter and for its properties 
thus becomes of crucial importance for building a fundamental theory of nature. Here we present 
a new universal relation, satisfied by matter distributions at all observed scales, and show its amaz- 
ingly good and detailed agreement with the predictions of the most up-to-date pure dark matter 
simulations of structure formation in the Universe [5-7]. This behaviour seems to be insensitive 
to the complicated feedback of ordinary matter on dark matter. Therefore, it potentially allows to 
compare theoretical predictions directly with observations, thus providing a new tool to constrain 
the properties of dark matter. This work extends the previous analysis [8-10] to a larger range of 
masses, demonstrates a different scaling law, and compares it with numerical simulations. Such 
a universal property, observed in structures of all sizes (from dwarf spheroidal galaxies to galaxy 
clusters), is difficult to explain without dark matter, thus providing new evidence for its existence. 

Current cosmological data, including the earliest stage of structure formation, are well-described by 
the "concordance model" (ACDM) [11]. The microscopic properties of the dark matter particles can 
strongly affect the formation of cosmic structures, originating from tiny preturbations in the primordial 
density field. Within the ACDM the dark matter is assumed to be non-interacting and to have no pri- 
mordial velocity - cold dark matter (CDM). To describe the formation and properties of virialized DM 
objects at a later, non-linear, stage of structure formation, numerical N-body simulations have been exten- 
sively used [12]. It was demonstrated that within the CDM model spherically averaged density profiles 
of DM halos of all scales are self-similar and universal [13, 14]. 

Many works have addressed the issue of measuring the DM distribution in observed objects (see 
e.g. [1-3, 15-17]). Two classes of profiles have been used to fit observational data: profiles with singular, 
cusp-like behaviour at small radii r (e.g. NFW profile [13]) or profiles that tend to a constant central 
density p c - cored profiles, such as e.g. pseudo-isothermal profile [18] (ISO) or Burkert profile [15] 
(BURK). Both profile families are described by two basic parameters: a characteristic scale r* at which 
the inner slope of the DM density profile changes towards its outer asymptotic and an average density p+ 
within this radius. 
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Unfortunately current data do not allow us to determine, in a conclusive way, the presence of cores or 
cusps in the observed density distributions. Often the same data set can be equally well fitted by profiles 
of different type [16, 17, 19-21]. The deviation from predictions of pure CDM models can occur for 
several reasons: for example, a significant influence of baryons; the microscopic properties of dark matter 
particles different from those of CDM (non-gravitational self-interaction or interaction with baryons, non- 
zero primordial velocities). Moreover baryons can contribute significantly (and may even dominate) the 
total mass profile in the inner regions, and thus pure CDM predictions are difficult to test. 

Recently the works [9, 10] analysed rotation curves and weak lensing data for a sample of dwarf, 
spiral and elliptical galaxies fitted by the Burkert profile and demonstrated that the mean dark matter 
surface density (S) = p*r* remains constant for all these galaxies. In this work we extend the analysis 
of [9, 10] to galaxies and galaxy clusters. We have complied a catalog of DM distribution in various 
celestial objects from more than 50 scientific publications. After applying uniform selection criteria to 
this catalog (see Method summary), we were left with 805 Dark Matter profiles (490 NFW, 285 ISO 
and 30 BURK) from 289 unique objects: 124 spiral galaxies, 1 1 dwarf spheroidals (dSphs), 10 elliptical 
galaxies, 25 galaxy groups and 121 galaxy clusters. To properly compare DM distributions, fitted by 
different density profiles, we introduce a dark matter column density, averaged over the central part of an 
object: 

5=4 f*rdr[dzp DM (V^T^)~^ , (1) 



Integral over z extends to the virial boundary of a DM halo. The definition ([]]) implies that S is pro- 
portional to the dark matter surface density within r* (S oc /9*r*). For distant objects S is defined via 
M cy i - mass within a cylinder of radius r*. Parameters of different profiles that fit the same DM den- 
sity distribution are related (for example, r s for NFW is equal to 6.1r c for ISO and equals to 1.6re for 
BURK). Choosing these values as r* in each case, one finds that the value of S for NFW (cuspy profile 
with r~ 3 asymptotic at large radii) and ISO (cored profile with r~ 2 behaviour at large scales) differ by 
less than 10% (the difference between NFW and BURK is ~ 2%), see Supporting Information for de- 
tails. Thus, the DM column density S is insensitive to the type of DM density profile, used to fit the same 
observational data. 

Our final dataset spans more than four orders of magnitude in r* (0.2 kpc < r* < 2.5 Mpc) and about 
eight orders of magnitude in halo massesQ The observational data can be fit by a single power-law: 

lg5 = °- 211g lb^ + L79 (2) 

(with S in Mq pc~ 2 ). One could try to interpret the data presented in Fig. 1 in spirit of [8-10] (i.e. 
S = const), although with a higher value of S 260M Q pc~ 2 (as was originally suggested in [22]). 
The apparent trend at higher masses could then be attributed to systematic errors and observational bias. 
Indeed, the results are based on observational data of different nature and different quality and we do 
not consider in this work observational errors. However, both Mh a i and S were computed, using 3 
density profiles on average, which should diminish the errors. The scaling relation (f2j) is supported by the 
data, spread over many orders in magnitude in halo mass, and the shift of some data points should not 
significantly affect the systematic trend. 

Next, we compare our data with outcomes from cosmological N-body simulations within the ACDM 
model [6, 7]. For each simulated halo we compute M^\ , fit the particle distribution to the NFW density 
profile and calculate S using formula ((T). The observational data together with results from ACDM 
numerical simulations [6] is plotted of the Fig. 1. The black dashed-dotted line on this Figure is the 
S — Mhaio relation obtained from N-body simulations [6], using the WMAP fifth year [11] cosmological 
parameters. It fairly well reproduces the fit ©. Moreover, the pink shaded region (showing the 3er scatter 



'We use M200 as halo mass M^\ . A proper definition of M200 can be found e.g. in [6] or in the Supplementary Information. 
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in the simulation data) contains most of the observational points within the halo mass range, probed by 
simulations. Therefore, the observed Aihaio — S scaling coincides with the relation between the parameters 
of DM density profiles observed in numerical simulations for long time [5,6,13] over more than five orders 
of magnitude in mass. 

Dwarf spheroidal satellites (dSphs) of the Milky Way (orange diamonds on the Fig. 1) do not follow 
the relation (0. Recently the Aquarius project has produced a statistically significant sample of well 
resolved density profiles for satellite halos [7], making it possible to determine their — relation. 
Satellites were found to be more concentrated than isolated halos and thus have a higher value of S 
at fixed Afhalo- Fig. 1 shows that the S — M^\ relation for satellite halos (gray dashed line) from the 
Aquarius simulation [7] reproduces well the data on dSphs. 

The fit to the data without the dSphs has the slope » 0.23, much better quality of fit, and coincides 
extremely well with the results of N-body simulations [6] for isolated halos (black dashed-dotted line on 
Fig. 1). At masses below 10 10 M Q no isolated halos were resolved in [6] and a simple toy model [6,23] 
was used to predict the relation between parameters of NFW profile in a given cosmological scenario. 
The model (dotted line in the Fig. 1) fits well the results for the few spiral galaxies in this range. Thus 
the agreement between observations and predictions from ACDM extend over more than eight orders of 
magnitude in mass. 

Comparison of our data with theoretical predictions (N-body simulations in our case) indicates that, 
despite the presence of various systematic errors in the data, the DM distributions in the observed objects 
exhibit a universal property - a systematic change of the average column density S as a function of the 
object mass (<S oc M^ a , relation (f5]i). This is different from the flat S = const dependence, previously 
suggested [8-10]. The latter is based, in our view, on a confusion between the properties of isolated and 
non-isolated halos. Excellent agreement with pure DM simulations suggest also that the observed scaling 
dependence is insensitive to the presence of baryons, details of local environment, formation history. 

The relation (0 can be used to search for deviations from CDM model (e.g. warm DM models [24]) 
or modifications of gravity at large scales [25]. This motivates dedicated astronomical observations with 
all the data processed in a uniform way. Studies of the galaxies with the masses below 1O 1O M0 and 
galaxy clusters would be especially important. 

Various scaling relations are known in astrophysics ("fundamental plane relation" for elliptical galax- 
ies [26], "Tully-Fisher relation" for spiral galaxies [27], etc.). The relation discussed in this Letter differs 
in one crucial aspect: it extends uniformly to all classes of objects at which DM is observed. It would be 
very difficult to explain such a relation within Modified Newtonian dynamics [28] theory considered as 
an alternative to DM. That is why this relation, further confirmed, studied and understood analytically, 
may serve as one more evidence of the existence of Dark Matter. 

If DM particle posses a radiative decay channel (as predicted by several particle physics models), a 
possible signal in X-ray or 7-ray observations would be proportional to the DM column density averaged 
within the instrument's field-of-view. It then follows from our analysis than decaying DM would produce 
a unique all-sky signal, with a known slow-varying angular distribution. Such a signal can be easily 
distinguished from any possible astrophysical background and therefore makes the astrophysical search 
for decaying DM another type of a direct detection experiment [22, 29, 30]. 

Method summary 

We collected from the literature 1095 DM density profiles for 357 unique objects ranging from dSphs to 
galaxy clusters. For each DM profile found in our sample we have applied uniform selection criteria: 

- If for an object several independently determined profiles were available and all of them but one agreed 
in the values of r+ and within a factor of 5, we rejected the outlier. 

- For some objects the best-fit value of the characteristic radius r* was extrapolated well outside the 
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Figure 1 : Column density S as a function of halo mass Mhaio- The black dashed-dotted line is the 
S — Mhaio relation obtained from N-body simulations [6], using the WMAP fifth year cosmological 
parameters [11]. The shaded region shows the 3cr scatter in the simulation data. The vertical lines 
indicate the mass range probed by simulations. The dotted line is the theoretical prediction from the toy 
model for isolated halos [6,23]. The gray dashed line showns the results from the Aquarius simulation 
for satellite halos [7]. 
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region covered by the observational data, i?data- In this case the parameters of the density profile had 
extremely large uncertainties. We have thus rejected objects with r* < 2.75i?data- 

- We rejected profiles for which the uncertainty in any quoted parameter (r* or p*) was higher than a 
factor of 10. 

- For objects with more than one profile selected, the average value of S and Mhaio wa s used in the 
subsequent analysis. 

- When processing the data of N-body simulations we used the fit of particle distribution by the NFW 
density profile and computed S, using equation ([l}. 

- If the same observational data is fit by several different DM profiles (e.g. NFW, ISO, and BURK), one 
can then find a relation between characteristic scales r* and densities p+ of these profiles. Provided 
such a relation holds, the difference between the column densities 5 NFW , 5 BU rk and S iS o turns out to 
be less than 10%. Qualitatively, this can be understood as follows: to explain the same velocity data, 
two DM profiles should have roughly the same mass within some radius Rq. If both profiles happen to 
have the same behaviour at large distances, their S values, averaged over Rq will be essentially equal. 
This explains the use of S as a characteristic of DM halos. 

Supplementary information is presented in the Appendix. 
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Supplementary Information 

In this Supplementary Information Section we demonstrate that the average DM column density S is 
independent on the choice of the particular DM density profile; provide more details about the sample 
that was used and its comparison with simulations. 

A Relation between parameters of DM density profiles 

In this work we concentrated on three popular choices for DM density profile p(r). Numerical (N-body) 
simulations of the cold DM model have shown that the DM distribution in all relaxed halos can be fitted 
with the universal Navarro-Frenk- White (NFW) profile [13]: 

( \ — Ps r s 

Pnfw(7) — —p. : -, y2 \?) 
r(l + r/r s y 
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parametrised by p s and r s . A more useful parametrization is in terms of the halo mass, M200, and the 
concentration parameter, c = i?2oo/ r s- Namely, -R200 is the radius at which the average DM density is 
200 times larger than the critical density of the universe p c rit- The halo mass M200 is the total mass of 
DM within this radius. The variables (p s ,r s ) and (M20Q, c) are thus connected as follows: 



/(c)Pcrit 

( 3M 200 

^OO^PcritC 3 



1/3 



(4) 



/(c) = - 



3 ln(c + l) -c/(c+l) 

The equations (|4|i allow to determine S oc r*p*, knowing halo mass and concentration parameter M20O) c. 

The Burkert profile [15] has been shown to be successful in explaining the kinematics of disk systems 
(e.g. [31]): 



Pburk ( r ) 



PBT% 



{tb + r)(r 2 B + r 2 ) ' 



(5) 



Another common parametrization of cored profiles is given by the pseudo-isothermal profile [18] 

Pc 



1 + r 2 /r 2 



(6) 



The quantity S(R) can be calculated analytically for all these choice of p(r). For example, for the 
pseudo-isothermal profile one obtains: 



<5iso (R) 

For the NFW density distribution (0: 



^p c r 2 
K 2 



VR 2 



R 2 



arctan y/R 2 /r 



1 



log 



R_ 

2r, 



(7) 



(8) 



Notice, that this expression is real for both R > r s and R < r s . The corresponding expression for the 
BURK is rather lengthy and not very illuminating. 



A.l Dependence of S on the inner slope of density profile. 

In order to equally well fit the same rotation curve data, two DM profiles should have roughly the same 
mass within some radius Rq, determined by the observational data. If both profiles happen to have the 
same behaviour at large distances, their S, averaged over Rq will be essentially equal (as it is determined 
by the sum of the masses inside the sphere Rq and in the outside of the cylinders, where the mass is 
dominated by the large r asymptotics)|3 In reality the situation is of course more complicated, one has to 
take into account the influence of baryons, the span of radii at which the data exists, etc. 

We conservatively estimate the difference of column densities between a cusped and a cored profile 
as follows. We take the NFW density profile © as a representative of the cusped profile and its "extreme 
cored" counterpart j o cor e (r) defined as follows: 

n (r\ - J ^pw(r), r > R ™ 

I PNFwi-floj, r S -n-0 
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Supplementary Figure 2: The ratio of average column densities of the extreme cored and NFW 
profiles as a function of R (c.f. Eq. ©)■ 



The column densities of these two profiles, averaged within Rq, differ only because the initial mass inside 
a sphere with radius Rq for the cored profile (O diminishes as compared to the NFW case. 

The resulting ratio of DM column densities is shown in the Supplementary Fig. [2] as a function of 
averaging radius Rq. In particular, for Rq = r s this ratio is 64%, for Ro — 2r s it equals to 53% and for 
Ro = 3r s it drops to 47%. This implies that the difference of DM column densities between the cusped 
(NFW) and the extreme version of the cored profile (|9]i is within 50% for realistic averaging radii Ro 
(usually Ro ~ 1 — 3r s ). This difference is small compared to the intrinsic scatter expected on a object 
by object basis and well below the observational uncertainties on the parameters describing the density 
profile. This makes S a very robust quantity to compare observed properties of DM halos and results 
from numerical simulations and, consequently, test the prediction of the CDM model. 

The rotation curve of a galaxy is often fitted by several DM profiles (e.g. ISO and NFW). Let us 
analytically establish the relation between parameters of several profiles, fitting the same rotation curve. 
To this we take an ISO density profile and generate according to it the circular velocity profile v\ (r), with 
r in the range r c < r < 15r c |j Then we fit these data using an NFW profile (see Fig.[3j left). We find the 
following relations between the parameters of the two profiles: 

NFW vs. ISO : r s ~6.1r c ; p„ ~ 0.11 p c . (10) 

The corresponding rotation curves and density profiles are shown in Supplementary Fig. [3] 

Let us now compare the column densities for NFW and ISO profiles, whose parameters are related 
via Eq. dlOb . Results as a function of radius R are shown on Fig. [4] In particular, one sees that for R = r s 



Siso(6r c ) 



One may be surprised that the cusped profile leads to the smaller column density than the cored one (as 
Eq. (fTTT i demonstrate). This result however, can be simply understood. We match the velocity profiles for 



2 We will see below (Eq. )13t ). that this is indeed the case for NFW and Burkert profiles. 
3 The final result is not sensitive to the exact choice of this range. 
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Supplementary Figure 3: Comparison of NFW (red) and ISO (green) profiles for the simulated ro- 
tation velocity points. Left panel: the velocity data (black points, in units of GnPc^) is generated, 
assuming the ISO profile and fitted with the NFW profile. The parameters of the corresponding NFW 
profile (in the units of r c , p c ) are given by eq. (1 1 Ob in the text. Right panel: comparison of the density 
profiles with parameters, related by (fTOb . The ir-axis is in the units of isothermal core radius r c . 



the NFW and ISO at some off-center distances i?o ~ 2r s , by demanding that the mass inside this sphere 
is the same for both profiles. The ISO profile is shallower in the outer regions than the NFW one. The 
ratio between the mass inside a sphere of the radius i?o and a cylinder with base radius i?o is equal to 0.58 
at i?o = 6r c for ISO profile, while it is 0.63 at Rq = r s for the NFW profile . Thus the mass in the outer 
part of a cylinder is larger for the shallower ISO profile than for the cuspy NFW one, which explains the 
result ( fTTT >. 

It is clear from previous considerations that 5 NHW and S BURK (similarly matched) should be essentially 
identical, as both profiles have identical behaviour at r — > oo. Indeed, in the case of the NFW and Burkert 
profiles the relation between their characteristic parameters is given by 

NFW vs. BURK : r s ~ 1.6r B ; p s ~ 0.37p B (12) 

which leads to 

<w^) ^ 098 (13) 



5burk (1.6t 



Finally, it should be noticed that we assume an infinite extension for Dark Matter halos, when com- 
puting the column density. However, the integrals in (O are convergent at large off-center distances and 
therefore the details of the truncation of the DM distributions for R > -R200 do not affect the value of S 
by more than 10%. 



B Data analysis 

We have collected from the literature 1095 DM profiles for 357 objects (from dwarf spheroidal galaxies 
to galaxy clusters, see Table[T]below). For each DM profile in our sample we have performed a number of 
checks. Those profiles have not passed these checks were rejected from subsequent analysis. As a result 
of the selection process we were left with 805 DM profiles for 289 objects. The results of the selection 
are shown on the Supplementary Figure |6]below. 

- When analysing the data, we realized that for some objects the value of r* lies well outside the region 
covered by the observational data, i?data- Such objects systematically show extremely high values of 
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Supplementary Figure 4: Comparison of the column densities of NFW and ISO profiles. Profiles 
describe the same data and their parameters are related via (llOl l. The column density is averaged within 
various radii R. Dashed vertical line marks R = r s = 6.1r c . 
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Supplementary Figure 5: Relation between parameters of NFW and ISO profiles in observed objects. 

For objects for which both NFW and ISO fits of velocity rotation curves were available, we plot the ratios 
r s /r c and p s j p c . The maximum of the histogram lies in a region ( 1 1 0b . 
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r*. For example, we found 37 galaxy profiles having r* > 100 kpc, while their kinematic data usually 
extends only up to i?data ~ 10 — 30 kpc. 

Therefore, we select only Dark Matter profiles having r* < 2.75i?d a ta- The coefficient 2.75 is justified 
by the following argument. The circular velocity in an NFW halo is given by 



v 2 c (r) 



GjvM nfw (< r) 



A 3 SV T, > r+r B 



(14) 



For r <C r s this function can be approximated as 



v c (r) w 2nG N p s r s [ r ~ — + 



4r 2 



(15) 



In the part of the velocity curve where is much less than the errors on the velocity dispersion one 
cannot reliably determine r s and p s (since v%(r) is indistinguishable from a straight line, proportional 
to p s r s ). It is important to have data points in the region where the contribution of the quadratic term 
becomes noticeable to reliably extract both NFW parameters. We chose to set 2.75i?data > ?V> which 
corresponds to a ~ 50% contribution from the second (quadratic) term to v%(r). Similar criteria are 
used for ISO and BURK profiles. This reduces the number of considered profiles from 1095 to 891. 

- For 76 objects both NFW and ISO (or BURK) Dark Matter profiles were available. For these objects 
we checked the relation between the parameters of these profiles against the results shown in Eq. ( fTOb 
(or Eq. (fT2l for BURK). Results of this comparison for the NFW and ISO profiles are shown in Fig. [5] 
This figure shows that there is indeed a maximum in the region defined by Eq. ( 1 1 0b but also that the 
scatter around this maximum is pretty large and that the difference between measured and expected 
ratios of NFW and ISO parameters can be as large as a factor of ten. Therefore we decided to exclude 
from our sample all objects with a ratio p s /p c , r s /r c , (or p s /ps, r s /rs for BURK profiles) larger than 
a factor 5 with respect to the theoretical prediction shown in Eq. ( llOl i or (fT2l . 

- Finally, in several cases parameters of DM density profiles were quoted with very large uncertainties. 
We decided to select only those profiles for which the ratio between the ler upper and lower bounds of 
quoted parameters (radius r* or the density p+) was smaller than a factor of 10. 

To compare the S — A/haio relation for selected objects with N-body simulations, we used the results 
from [6]. This suit of ACDM numerical simulations probed the halo mass range 10 10 — 10 15 M Q . For 
each simulated halo of [6] we computed Mhaio> fit the particle distribution to the NFW density profile 
and calculate S using Eqs.(0]i and the definition (Q]). The observational data together with results from 
simulations is plotted of the Supplementary Figure[7] The small scatter of the simulation points at Mhaio ^ 
10 14 M Q is explained by the finite size of the simulation box. The simulations with the large box size 
(e.g. [5]) verify that the scatter does not reduce at large masses (c.f. the pink shaded region on the 
Figure 1). 

Supplementary Table 2: Halo mass M200 and DM column density S for galaxy clusters. 
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2A 0335+096 


1.44e+15 




2.67 




Abell 1068 


5.64e+14 




2.59 




Abell 1 19 


2.64e+15 




2.74 




Abell 133 


4.35e+14 




2.66 




Abell 1367 


6.89e+14 




2.72 




Abell 1413 


9.87e+14 




2.78 




Abell 160 


9.73e+13 




2.63 




Abell 1656 (Coma) 


1.72e+15 




3.01 




Abell 168 


3.19e+14 




2.64 




Abell 1689 


1.56e+15 




3.28 




Abell 1703 


1.25e+15 




3.28 




Abell 1763 


9.92e+14 




2.51 




Clusters (continued on next page) 
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Object name 


M 20 o[M e ] 
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_ M pc- 2 _ 




_ M(T) pc~ 2 _ 




Abell 1795 


1.76e+15 


2.79 


Abell 1835 


1.27e+15 


2.68 


Abell 194 


8.90e+13 


2.53 


Abell 1983 


1.57e+14 


2.39 


Abell 1991 


1.64e+14 


2.69 


Abell 2029 


3.45e+15 


2.90 


Abell 209 


5.26e+14 


2.47 


Abell 2142 


5.98e+15 


2.65 


Abell 2199 


9.57e+14 


2.81 


Abell 2204 


7.24e+14 


2.90 


Abell 2218 


7.88e+14 


2.98 


Abell 2219 


1.65e+15 


2.77 


Abell 2256 


3.84e+15 


2.58 


Abell 2319 


4.75e+15 


2.68 


Abell 2390 


1.54e+15 


2.82 


Abell 2537 


8.18e+14 


2.82 


Abell 2597 


2.94e+14 


2.73 


Abell 262 


1.89e+14 


2.68 


Abell 2667 


1.37e+15 


2.64 


Abell 267 


1.74e+14 


2.51 


Abell 2717 


1.66e+14 


2.44 


Abell 3112 


2.84e+14 


2.83 


Abell 3266 


5.17e+15 


2.62 


Abell 3376 


2.13e+14 


2.75 


Abell 3526 


9.84e+14 


2.87 


Abell 3562 


9.33e+14 


2.86 


Abell 3571 


5.17e+15 


2.78 


Abell 3627 


1.43e+15 


2.77 


Abell 370 


1.25e+15 


2.88 


Abell 383 


6.19e+14 


2.56 


Abell 426 


1.66e+15 


2.97 


Abell 478 


1.42e+15 


2.76 


Abell 496 


1.64e+15 


2.99 


Abell 511 


9.21e+14 


2.90 


Abell 520 


1.29e+16 


3.07 


Abell 539 


2.83e+14 


3.22 


Abell 576 


6.47e+14 


3.12 


Abell 611 


7.43e+14 


2.82 


Abell 644 


6.96e+14 


2.72 


Abell 68 


3.71e+14 


2.54 


Abell 754 


5.32e+15 


2.62 


Abell 780 


9.96e+13 


3.02 


Abell 85 


2.38e+15 


2.74 


Abell 851 


9.85e+13 


2.42 


Abell 907 


6.24e+14 


2.79 


Abell 963 


6.78e+14 


2.87 


Abell SI 101 


1.76e+14 


2.79 


Abell S400 


3.20e+15 


3.04 


Abell S540 


1.69e+14 


2.71 


AC 114 


1.27e+15 


2.67 


AWM 4 


1.50e+14 


2.71 


CL 0016+1609 


1.31e+14 


2.19 


CL 2244-0221 


4.30e+14 


2.65 


CL J1226.9+3332 


7.04e+14 


3.29 


B2002a 15 


2.87e+14 


2.65 


MACS J001 1.7-1523 


1.14e+15 


2.71 


MACS J0159.8-0849 


1.25e+15 


2.94 


MACS J0242.6-2132 


5.28e+14 


2.98 


MACS J0329.7-0212 


7.19e+14 


2.86 


MACS J0429.6-0253 


4.1 le+14 


3.04 


MACS J0744.9+3927 


1.01e+15 


2.93 


MACS J0947. 2+7623 


1.16e+15 


2.97 


MACS Jl 115.8+0129 


4.32e+15 


2.59 


MACS J131 1.0-03 11 


6.72e+14 


2.82 


MACS J1411. 3+5212 


7.16e+14 


3.13 


MACS J1423. 8+2404 


5.79e+14 


3.13 


MACS J1427. 6-2521 


3.07e+14 


2.95 


MACS J1532.9+3021 


9.40e+14 


2.86 


MACS J1621. 6+3810 


7.58e+14 


2.99 


MACS J1720.3+3536 


9.67e+14 


2.85 


MACS J1931. 8-2635 


1.67e+15 


2.78 


MACS J2229.8-2756 


3.10e+14 


2.97 


MKW 3S 


1.10e+14 


2.34 


MKW 9 


1.19e+14 


2.54 


MS 0015.9+1609 


9.50e+15 


3.21 


MS 0116.3-0115 


1.18e+14 


2.47 


MS 0302.7+1658 


8.45e+14 


3.12 


MS 0451.6-0305 


1.68e+15 


3.03 


MS 0735.6+7421 


2.17e+15 


3.15 


MS 0839.9+2938 


5.99e+14 


2.93 


MS 1006.0+1202 


3.10e+15 


2.93 


MS 1008.1-1224 


3.39e+15 


2.99 


MS 1137.5+6625 


6.21e+14 


2.69 


MS 1224.7+2007 


9.01e+14 


3.37 


MS 1358.4+6245 


2.58e+15 


3.12 


MS 1455.0+2232 


1.37e+15 


3.38 


MS 1512.4+3647 


7.11e+14 


3.12 


MS 1910.5+6736 


8.59e+14 


2.80 


MS 2137.3-2353 


5.11e+14 


3.13 


PKS 0745-19 


2.14e+15 


2.96 


PKS 0745-191 


2.87e+15 


2.82 


RXC J0137.2-0911 


9.83e+13 


2.61 


RXC J0454.8-1806 


1.21e+14 


2.58 


RXC J1504.1-0248 


1.95e+15 


2.84 


RX J0056.9-2740 


6.45e+13 


1.95 


RX J0419.6+0225 


3.06e+13 


2.87 


RX J0439.0+0520 


4.39e+14 


2.94 


RX J0439.0+0521 


3.85e+14 


2.86 


RXJ1023.6+0411 


2.88e+15 


2.62 


RX J1223. 0+1037 


3.31e+13 


2.62 


RX J1347.5-1 145 


2.36e+15 


3.24 


RX J1416.4+2315 


3.00e+14 


3.14 


RX J1504.1-0248 


1.80e+15 


2.77 


RX J2129.6+0005 


7.15e+14 


2.71 


SDSS J1004+4112 


3.02e+14 


2.87 


Triangulum 


3.32e+15 


2.91 


Virgo 


4.16e+14 


2.35 


ZWCL 0024.0+1652 


3.49e+14 


2.80 


ZWCL 1321.4+1358 


1.50e+13 


2.68 








Clusters (end) 
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Supplementary Figure 6: Selected DM density profiles. DM column density S as a function of the 
characteristic radius r* for 1095 density profiles. Filled (open) shapes denote accepted (rejected) density 
profiles, according to our selection criteria. Notice, that unlike Figure 1, an object may appear several 
times in this Figure if more than one profile was available. 



Type 


Number of objects 


References 


List of selected objects 


Galaxy clusters 


130 


[32^14] 


Table |2] 


Galaxy groups 


26 


[45-47] 


Table [3] 


Elliptical galaxies 


10 


[48-52] 


Tableg] 


Spiral galaxies 


180 


[18,31,53-79] 


Tabled 


Dwarf spheroidal galaxies 


11 


[16,79-83] 


Tableg] 


total 


357 







Supplementary Table 1: Observational data. The table lists the types of objects; references used to 
collect the observational data; and the final list of selected objects. 
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Supplementary Figure 7: DM Column density as a function of the halo mass. Similar to the Figure 1, 
we plot 289 objects, selected in Section [B] above (coloured shapes) superimposed on the simulation data 
for isolated halos [6] (open black circles). 
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Supplementary Table 3: Halo mass M200 and DM column density S for galaxy groups. 



Object name 


M200 [M Q ] 


lo SlO 


5 




Object name 


M 20 o[Mq] 


lo Sl0 


S 




_ M(7) pc~ 2 _ 


_ M R pc~ 2 _ 


Abell 1275 


8.20e+13 


2.57 


Abell 1692 


1.09e+14 


2.77 


Abell 2462 


l.lle+14 


2.51 


Abell 2550 


8.93e+13 


2.64 


Abell 3880 


2.03e+14 


2.69 


ESO 306-G017 (group) 


1.74e+14 


2.73 


ESO 351-G021 (group) 


3.66e+13 


2.18 


ESO 552-G020 (group) 


1.31e+14 


2.60 


IC 1860 (group) 


5.02e+13 


2.58 


MKW4 


5.97e+13 


2.68 


NGC 1407 (group) 


6.13e+13 


2.52 


NGC 1550 (group) 


3.61e+13 


2.73 


NGC 2563 (group) 


2.51e+13 


2.51 


NGC 4104 (group) 


6.18e+13 


2.55 


NGC 4325 (group) 


3.57e+13 


2.64 


NGC 5044 (group) 


3.67e+13 


2.64 


NGC 5098 (group) 


2.52e+13 


2.42 


NGC 5129 (group) 


1.51e+13 


2.69 


NGC 533 (group) 


2.31e+13 


2.84 


RGH80 


2.80e+13 


2.53 


RXJ1022.1+3830 


1.00e+14 


2.58 


RXJ1158.1+5521 


1.05e+14 


2.43 


RXJ1159.8+5531 


8.33e+13 


2.75 


UGC 5088 (group) 


1.93e+13 


2.36 


UGC 842 (group) 


7.00e+13 


2.77 








Galaxy Groups 



Supplementary Table 4: Halo mass M200 and DM column density S for elliptical galaxies. 



Object name 


M 20 o[M e ] 


l°gl<) 


S 




Object name 


M 2 oo[M Q ] 


lo §10 


S 




_ M fl pc~ 2 _ 




_ M fl pc~ 2 _ 




M49 (NGC 4472) 


2.85e+13 




2.67 


M60 (NGC 4649) 


3.05e+13 




2.99 


NGC 4125 


5.32e+12 




2.26 


NGC 4261 


4.77e+13 




2.01 


NGC 6482 


6.12e+12 




2.64 


NGC 720 


5.43e+12 




2.63 


Elliptical galaxies 



Supplementary Table 5: Halo mass M200 and DM column density <S for dwarf spheroidal 
galaxies (dSphs). 



Object name 


M 20 o[M©] 


!ogio 


s 




Object name 


M 2OO [M ] 


logio 


s 




_ M fl pc~ 2 _ 




_ Mq pc~ 2 _ 




Carina 


2.63e+()9 




1.64 


Coma Berenices 


1.00e+08 




1.85 


Draco 


1.25e+09 




1.97 


Fornax 


3.04e+09 




1.54 


Leo II 


7.17e+08 




1.69 


Ursa Major II 


1.27e+08 




1.83 


Ursa Minor 


1.85e+09 




2.01 


Willman I 


7.61e+07 




1.98 


Dwarf spheriodal galaxies 
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Supplementary Table 6: Halo mass M200 and DM column density S for spiral, low-surface 
brightness (LSB) and dwarf galaxies (without dSphs). 



Object name 


M 20 o[M e ] 


lo SlO 


s 




Object name 


M 2O o[M ] 




s 












DDO 154 


5.52e+10 


1.44 


DDO 170 


6.82e+10 


1.60 


DDO 189 


4.26e+10 


1.85 


DDO 47 


2.13e+ll 


2.03 


DDO 52 


2.00e+10 


1.89 


DDO 64 


5.28e+10 


1.89 


ESO 116-G12 


5.14e+10 


1.65 


ESO 287-G13 


1.36e+09 


0.95 


ESO 79-G14 


4.41e+08 


1.08 


IC 1401 


1.54e+12 


1.31 


IC 2574 


2.28e+ll 


1.59 


KK 98250 


4.57e+10 


1.70 


LSBCF563-1 


2.78e+ll 


1.95 


LSBC F563-V2 


2.41e+ll 


2.10 


LSBCF568-1 


4.35e+ll 


2.28 


LSBC F568-3 


4.66e+ll 


1.89 


LSBC F568-vl 


2.37e+ll 


2.25 


LSBC F57 1-8 


5.74e+ll 


2.25 


LSBCF574-1 


2.15e+ll 


2.05 


LSBCF583-1 


2.56e+ll 


1.75 


LSBC F583-4 


9.37e+10 


1.82 


M31 


8.30e+ll 


2.63 


M33 (NGC 598) 


6.48e+ll 


1.93 


Milky Way 


1.03e+12 


2.40 


NGC 100 


2.50e+ll 


2.04 


NGC 1003 


5.30e+ll 


1.48 


NGC 1085 


4.37e+12 


2.35 


NGC 1087 


7.06e+ll 


2.10 


NGC 1090 


3.01e+08 


1.06 


NGC 1353 


1.27e+12 


2.34 


NGC 1421 


1.60e+12 


1.93 


NGC 1560 


7.91e+ll 


1.67 


NGC 224 


1.18e+13 


1.97 


NGC 2366 


2.51e+10 


1.79 


NGC 2403 


3.80e+ll 


2.11 


NGC 247 


4.47e+12 


1.88 


NGC 2532 


4.39e+ll 


2.18 


NGC 2537 


7.85e+10 


2.29 


NGC 2552 


1.44e+12 


2.38 


NGC 2608 


3.10e+ll 


1.94 


NGC 2708 


1.86e+12 


2.63 


NGC 2841 


1.87e+12 


2.66 


NGC 2903 


6.96e+ll 


2.43 


NGC 2977 


4.54e+ll 


2.47 


NGC 300 


8.95e+ll 


1.87 


NGC 3026 


1.39e+12 


2.43 


NGC 3031 


2.52e+ll 


2.70 


NGC 3109 


1.16e+ll 


1.30 


NGC 3198 


3.64e+ll 


2.04 


NGC 3200 


2.16e+12 


2.53 


NGC 3274 


5.50e+10 


2.46 


NGC 3319 


8.11e+ll 


1.74 


NGC 3346 


4.04e+ll 


2.50 


NGC 3495 


1.68e+13 


3.33 


NGC 3521 


5.80e+ll 


2.15 


NGC 3621 


6.05e+ll 


1.93 


NGC 3741 


1.99e+10 


1.73 


NGC 3893 


1.98e+12 


2.84 


NGC 3898 


9.42e+12 


2.48 


NGC 4190 


5.20e+ll 


1.59 


NGC 4236 


1.95e+ll 


1.73 


NGC 4242 


1.08e+12 


1.80 


NGC 4258 


3.04e+12 


1.92 


NGC 4288 


1.00e+ll 


2.26 


NGC 4395 


6.40e+10 


2.06 


NGC 4455 


1.00e+ll 


1.78 


NGC 4494 


8.65e+ll 


1.87 


NGC 45 


8.25e+ll 


1.97 


NGC 4605 


1.06e+ll 


2.33 


NGC 4635 


5.21e+ll 


2.53 


NGC 4736 


2.30e+ll 


2.48 


NGC 4800 


3.51e+ll 


2.70 


NGC 5023 


1.35e+ll 


2.14 


NGC 5033 


1.25e+12 


1.93 


NGC 5055 


2.17e+12 


1.86 


NGC 5194 


1.56e+12 


2.82 


NGC 5204 


9.53e+10 


2.15 


NGC 55 


1.00e+12 


1.73 


NGC 5585 


3.40e+ll 


1.85 


NGC 5608 


4.97e+10 


1.68 


NGC 5727 


5.41e+12 


2.06 


NGC 5949 


1.72e+ll 


2.33 


NGC 5963 


6.08e+ll 


1.76 


NGC 6015 


1.25e+12 


2.49 


NGC 6503 


6.90e+ll 


1.98 


NGC 6682 


5.22e+l() 


1.78 


NGC 6689 


2.76e+ll 


2.12 


NGC 6946 


6.19e+12 


2.06 


NGC 7137 


4.45e+10 


2.12 


NGC 7217 


2.47e+12 


2.69 


NGC 7339 


3.54e+ll 


1.89 


NGC 753 


1.92e+12 


1.82 


NGC 7541 


2.52e+12 


2.17 


NGC 7664 


1.26e+12 


2.03 


NGC 7793 


3.68e+ll 


2.16 


NGC 925 


6.16e+ll 


1.81 


NGC 959 


4.87e+10 


2.60 


UGC 02259 


9.49e+ll 


1.76 


UGC 10310 


1.81e+ll 


1.71 


UGC 11557 


4.47e+ll 


1.84 


UGC 11707 


2.11e+ll 


2.03 


UGC 11820 


1.99e+ll 


2.43 


UGC 11861 


8.14e+ll 


2.21 


UGC 12060 


5.37e+ll 


1.98 


UGC 1230 


1.47e+ll 


2.01 


UGC 12632 


4.55e+10 


1.92 


UGC 12732 


1.47e+ll 


1.81 


UGC 128 


3.84e+ll 


2.08 


UGC 1281 


9.96e+10 


1.74 


UGC 12810 


1.78e+12 


2.40 


UGC 1551 


4.98e+10 


1.82 


UGC 2259 


5.37e+10 


2.35 




Spiral, LSB and dwarf 


galaxies (without dSphs) (continued on next page) 
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Object name 


M200 [Mq] 


lo gl0 


5 




Object name 


M 20 o[M e ] 


lo SlO 


S 




_M R pc- 2 _ 


.Mr pc- 2 _ 


UGC 2885 


8.22e+12 




1.86 


UGC 3137 


2.10e+ll 




2.02 


UGC 3371 


1.81e+ll 




1.75 


UGC 4173 


4.62e+10 




1.23 


UGC 4499 


1.03e+ll 




1.85 


UGC 477 


4.00e+ll 




1.98 


UGC 5005 


4.52e+ll 




1.54 


UGC 5750 


2.50e+ll 




1.46 


UGC 628 


2.73e+ll 




2.23 


UGC 6446 


5.83e+10 




1.84 


UGC 731 


4.42e+10 




2.13 


UGC 7699 


1.44e+ll 




1.70 


UGC 9211 


2.84e+10 




1.99 










Spiral, LSB and dwarf galaxies (without dSphs) (end) 
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